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Abstract 

Multiple emission of intermediate-mass fragments has been studied for the collisions 
of p, 4 He and 12 C on Au with the 4vr setup FASA. In the case of 12 C(22.4 GeV)+Au 
and 4 He(14.6 GeV)+Au collisions, the deviations from a pure thermal break-up are 
seen in the energy spectra of the emitted fragments: the spectra are harder than 
calculated and than measured in p-induced collisions. This difference is attributed to 
a collective flow with the expansion velocity on the surface about 0.1 c (for 12 C+Au 
collisions). 
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Nuclear multifragmentation is a new decay mode of excited nuclei character- 
ized by the emission of Intermediate Mass Fragments (IMF, 3 < Z < 20). The 
development of this field has been strongly stimulated by the idea that this 
process might be related to a liquid-gas phase transition in nuclear matter. A 
recent status on multifragmentation can be found in Ref. [1]. 

Highly excited nuclei can be formed by heavy ion collisions at intermedi- 
ate energies, but there are advantages in producing the excitation by light 
relativistic projectiles. In the first case, nuclear heating is accompanied by 
compression, fast rotation and shape distortion which may cause dynamic ef- 
fects in the mult i- fragment disintegration and it is not easy to disentangle 
all these effects and extract information on thermodynamic properties of hot 
nuclear systems. The situation becomes more transparent if light relativistic 
projectiles are used. In this case, dynamic effects are expected to be negligible. 
Another advantage is that all the fragments are emitted by a single source: a 
slowly moving target remainder. Its excitation energy might be almost entirely 
thermal. Light relativistic projectiles provide therefore an unique possibility 
to study "thermal multifragmentation". It has been shown that thermal mul- 
tifragmentation indeed takes place in collisions of light relativistic projectiles 
(p, 3 He, a and recently it and p) with a heavy target [2-12]. 

In this Letter we concentrate on studying energy spectra of the emitted frag- 
ments as they reflect, due to the "Coulomb law", the geometry and dynamics 
(expansion) of the emitting source. By comparing our results from p+Au [4] 
collision with those from reactions induced by a particles and 12 C projectiles 
with incident energies of (1 - 4) A-GeV, we evidence a transition from a pure 
statistical process ("thermal multifragmentation") to a behaviour reflecting 
dynamics. It will be shown that from the observed additional collective en- 
ergy a spatial distribution of the fragments at freeze out can be infered. 

The experiments were performed with beams from the JINR synchrophasotron 
in Dubna using the modified [13] 47r-setup FASA [14]. The device consists of 
two main parts : (i) Five AE (ionization chamber) x E (Si)-telescopes, which 
serve as a trigger for the read-out of the system allowing measurement of the 
charge and energy distributions of IMF's at various angles from 24° to 156° 
covering together a solid angle of 0.03 sr. (ii) The fragment multiplicity detec- 
tor consisting of 64 CsI(Tl) counters (with thicknesses around 30 mg-cm -2 ) 
which covers 89% of 4ir. This device gives the number of IMF's in the event and 
their spatial distribution. A self-supporting Au target of 1.5 mg/cm 2 thickness 
was located in the center of the FASA vacuum chamber (~ 1 m in diameter). 
The following beams were used: protons at energies of 2.16, 3.6 and 8.1 GeV, 
4 He at energies of 4 and 14.6 GeV and 12 C at 22.4 GeV. The average beam 
intensity was 7 ■ 10 s particles/spill for protons and helium and 1 • 10 8 parti- 
cles/spill for carbon projectiles with a spill length at 300 ms and a spill period 
of 10 s. See also [4] reporting on the p+Au experiment. 
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The kinetic energy of fragments is determined by four terms: thermal motion, 
Coulomb repulsion, rotation, and collective expansion energies of the system 
at freeze out, E = E th + E c + E rot + Ef low . The additivity of the first three 
terms is quite obvious. For the last term, its independence of others may be 
considered only approximately when the evolution of the system after the 
freeze-out point is driven only by the Coulomb force. The Coulomb term is 
significantly larger than the thermal one as was shown in Ref. [7] for 4 He 
(14.6 GeV)+Au collisions, where the Coulomb part of the mean energy of the 
carbon fragment is three times larger than thermal one using volume emission 
of fragments from a diluted system. 

The contribution of the collective flow for the p+Au collisions at 8.1 GeV 
incident energy was estimated to Vfi ow < 0.02 c in Ref. [4]. This was done 
by comparing the measured IMF spectra with the ones calculated within the 
framework of the Statistical Multifragmentation Model (SMM) [15] which in- 
cludes no flow. For heavy ion collisions, collective flow has been observed and 
it is the most pronounced in central Au+Au collisions [16]. In this respect, it 
is interesting to analyse energy spectra of fragments from He+Au and C+Au 
collisions looking for a possible onset of collective flow phenomenon. 

A comparison of the energy spectra of outgoing carbon isotopes from proton-, 
helium- and carbon-induced collisions on a Au target is presented in Fig. 1. 
The spectral shapes show an increase in the number of high-energy carbon 
fragments with the projectile mass. 

The reaction mechanism for light relativistic projectiles is usually divided into 
two stages. The first one is a fast energy-depositing stage, during which very 
energetic light particles are emitted and a nuclear remnant (target spectator) is 
excited. We use a refined version of the intranuclear cascade model (INC) [17] 
for describing the first stage. The second stage is described by the Statistical 
Multifragmentation Model (SMM), which considers multibody decay (volume 
emission) of a hot expanded nucleus. But such a two-stage approach fails to 
describe the observed IMF multiplicities as shown in Table 1. An expansion 
stage is inserted between the two parts of calculation. The excitation energies 
and the residual masses are then fine tuned [4] to get an agreement with the 
measured IMF multiplicities, i.e. the values for the residual (after INC) masses 
Ar and their excitation energies Er are scaled on an event-by-event basis. 
The average masses of nuclei which decay by multifragmentation is labelled 
A MF having mean excitation energies E MF as given in Table 1 together with 
the values of < M IMF >. The lines in Fig. 1 give the spectra calculated in 
the framework of this combined model INC+Expansion+SMM. The fragment 
energies are obtained by the multibody Coulomb trajectory calculations on an 
event-by-event basis. In the initial state all the charged particles are assumed 
to have a thermal velocity only (no flow). 
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The calculated carbon spectrum for p+Au collisions (at 8.1 GeV) is consistent 
with the measured one. A similar situation occurs with 4 He+Au collisions 
at 4 GeV, but not with 4 He(14.6 GeV)+Au and 12 C+Au interactions: the 
measured spectra are harder than the calculated ones. A careful inspection of 
the reaction induced with a particles shows that the maximum of the energy 
spectra obtained at 4 GeV is higher than obtained at 14.6 GeV. This is mainly 
caused by the higher Z values of the residual nuclei at the lower incident 
energy and is reproduced by the model calculations. The slope of the 14.6 
GeV data is much harder than for the lower incicent energy. This higher slope 
is underestimated by the model and this failure is most pronounced for C + 
Au collisions. An attempt to describe these spectra by a higher Coulomb field, 
caused e.g. by a more compact system at break-up fails. In such an approach 
the maximum shifts keeping the drop towards higher kinetic energies while 
the measured spectra show an increase at the higher kinetic energies without 
moving the maximum significantly. 

The trends from Fig. 1, i.e. increasing mean energies with increasing mass and 
increasing energy of the projectile, is seen for many emitted fragments. This 
observation is summarized in Fig. 2, which shows the mean kinetic energies 
per fragment nucleon as a function of the charge of the detected fragments Z. 
This figure shows a remarkable enhancement in the kinetic energies for light 
fragments emitted in He+Au and C+Au collisions as compared to the p(8.1 
GeV)+Au case. The calculated values of the mean fragment energies (shown 
by lines) are obtained with INC+Expansion+SMM model. The measured en- 
ergies are close to the calculated ones for p+Au collisions in the range of the 
fragment charges between 4 and 9. The experimental values for 4 He+Au and 
12 C+Au interactions, however, exceed the calculated ones, which are similar 
for all three cases. 

The observed deviation cannot be attributed to an angular momentum effect. 
The rotational energy E rot of a fragment with mass Ajmf can be estimated 
from the total rotational energy El of a system with mass number Ar using 
classical rotation: 

< E rot > /Ajmf = - < -± > (1) 
6 A R K sys 

where Rz is the radial coordinate of the IMF and R sys is the radius of the 
system. According to the INC calculations for C+Au collisions, the mean 
angular momentum L of the target spectator is ~ 36ft. It might be reduced 
by a factor of 1.5 due to the mass loss during expansion along the way to 
the freeze-out point. Finally < El > is estimated to be only 5 MeV and 
< E rot > /Aj MF ~ 0.04 MeV/nucleon, which is by an order of magnitude 
smaller than the energy enhancement for light fragments. We suggest that 
the observed enhancement is caused by the expansion of the system, which 
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seems to be radial, as a v ||-versus--u j_ plot (this will be subject of a forthcoming 
publication) does not show any significant deviation from circular symmetry. 

An estimate of the fragment flow energy may be obtained from the difference 
between the measured IMF energies < E > and those calculated (without any 
flow) < E > no fiow This difference for C+Au collisions is shown in the lower 
part of Fig. 2. A drastic decrease of the flow energy with increasing charge of 
the fragment is seen. 

In an attempt to describe the data, we modified the SMM code by including 
a radial velocity boost for each particle at freeze out, i.e. a radial expansion 
velocity was superimposed on the thermal motion in the calculation of the 
multibody Coulomb trajectories on an event-by-event basis. A self-similar ra- 
dial expansion of a spherical nucleus is assumed, where the local flow velocity 
is linearly dependent on the distance of the particle from the centre of mass. 
The expansion velocity of particle Z located at radius Rz is given by the 
following expression 

— # 

v now (Z) = v° flow ■ — — (2) 



where Vf low is the radial velocity on the surface of the system. The use of the 
linear profile of the radial velocity is motivated by the hydrodynamic models 
for an expanding hot nuclear system (see for example Ref. [18]). The value of 
vj low has been adjusted to 0.1 c in order to describe the mean kinetic energy 
measured for the carbon fragments. The results are presented in the lower 
part of Fig. 2 as dashed line calculated from the difference of the fragment 
energies obtained for v j iow =0.1 c and v® low =0. The data deviate significantly 
from the calculated values for Li and Be. This may be caused in part by the 
contribution of particle emission during the early stage of expansion from the 
hotter and denser system. It is supported by the fact that the extra energy of 
Li fragments with respect to the calculated value is clearly seen in Fig. 2 even 
for the proton-induced fragmentation, where no significant flow is expected. 
This peculiarity of light fragments has been noted already by the ISIS group 
for 3 He+Au collisions at 4.8 GeV [10]. 

For fragments heavier than carbon, the calculated curve in the lower part 
of Fig. 2 is higher than the data and only slightly decreasing with increasing 
fragment charge. The trend of the calculation is to be expected. The mean frag- 
ment flow energy is proportional to < R% > and this value changes only little 
with fragment charge in the SMM code due to the assumed equal probability of 
IMF distribution inside the available break-up volume. The difference between 
data and calculations shown in Fig. 2 indicates that a uniform density distri- 
bution is not fullfilled. A preference of finding heavier fragments more in the 
center would reduce their Coulomb energies shown in the upper part of Fig. 2, 
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and would increase the extracted quantity < E rneasured > — < E > no fiow 
leading then to a consistent description of a flow scenario with a non-uniform 
fragment distribution. The dense interior of the expanded nucleus may be 
favored for the appearance of larger IMF's, if fragments are formed via the 
density fluctuations. Indications for such an effect could already be drawn 
from the analysis of the mean IMF energies performed in Ref. [4] for proton 
induced fragmentation: the measured energies are below the theoretical curve 
for fragments heavier than Ne. This means that our procedure of estimating 
the flow energy is only reliable for fragments lighter than Ne where the model 
calculations fit the data of the fragment energies from p+Au collisions (no 
flow). For heavier fragments the flow energies are underestimated and should 
be considered as lower limits. The interesting feature of a reduction of the 
flow energy for heavier fragments is observed also for the central heavy ion 
collisions [19]. This effect is increasingly important at energies <100 A-MeV, 
and that is in accordance with our suggestion on its relation to the density 
profile of the hot system at freeze out. 

For the estimation of the mean fragment flow velocities < Vfi ow > the differ- 
ence between the measured IMF energies and calculated ones (without flow) 
has been used. The results are presented in Fig. 3. The values for Li and 
Be are considered as upper limits because of the possible contribution of the 
preequilibrium emission. The corresponding values of the fragment positions 
< Rz > I R sys for the freeze-out point obtained under assumptions of the lin- 
ear radial profile for the expansion velocity are given on the right-hand scale 
of this figure. The dotted line shows the mean radial coordinates of fragments 
according to the SMM code. As it has been noted above, the calculated val- 
ues of < Rz > /R sys are only slightly decreasing with Z as expected from a 
uniform density distribution, but in clear contrast to the data. Calculations 
with the MMMC model [20] give the same trend as the SMM model [21]. 

Effects of the radial collective energy for 1 A-GeV Au+C collisions (in inverse 
kinematics) were considered in [22] by analysing the transverse kinetic ener- 
gies. The mean radial flow velocities were estimated, but it had been done 
only for fragments with Z < 7. In this charge range our analysis gives slightly 
lower values of the mean expansion velocities as compared to Ref. [22]. Their 
interpretation of a time sequence of the emission acccording to the Z value is 
largely equivalent to our interpretation. 

The total expansion energy can be estimated by integrating the nucleon flow 
energy (taken according to Eq. (2)) over the available volume at freeze out. 
For an uniform system one obtains 

EfL = ^A R ■ m N (v° flow ) 2 (1 - r N /R sys f (3) 



6 



with the mass itin and the radius of a nucleoli. For 12 C+Au collisions it 
gives Ejf^ ~ 100 MeV, corresponding to the flow velocity at the surface of 
0.1c. Similar results are obtained for 4 He(14.6 GeV)+Au collisions. The total 
flow energy of the fragmenting nucleus is four times less than the thermal one 
estimated in the INC+Expansion+SMM approach. 

Concluding, the energy spectra of IMF's for reactions p+Au at 2.1, 3.6 and 8.1 
GeV, 4 He+Au at 4 and 14.6 GeV, 12 C+Au at 22.4 GeV are compared. While 
the fragment kinetic energies are well described within the SMM code for 
p+Au collisions, the model underestimates the kinetic energies of fragments 
from collisions induced by 4 He (14.6 GeV) and 12 C (22.4 GeV) projectiles. 
The additional energy is attributed to collective expansion. However, a linear 
flow profile fails to describe the variation of flow energies extracted from the 
measured spectra with the fragment charge. This discrepancy might be caused 
by the fact that the model assumes a uniform density distribution and, hence, 
a rather uniform probability distribution of forming fragments in the freeze- 
out volume. The data indicate that heavy fragments are preferentially located 
more in the interior of the nucleus. 

The presented study shows that in spite of the success of statistical multi- 
fragmentation models in describing the partitions, the freeze-out condition 
might be still too simplified. The energy spectra provide sensitive probes for 
the source configuration and the emission dynamics. The range of projectiles, 
from protons to light nuclei, seems to be quite attractive in this respect for 
showing a transition from "thermal break-up" to a disintegration possibly 
caused by rapide expansion likely together with a non-uniform density profile 
of the excited nucleus. 
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Proj 


Exper. 


Calculations 


Model 


(GeV) 




MjmF 


MjmF 


Ar 


Amf 


Er 


Emf 




2.16 


P 


1.7±0.2 


1.82 


189 


185 


310 


589 


INC+SMM 








1.69 


188 


183 


288 


564 


with expansion 


3.6 


P 


1.9±0.2 


2.52 


187 


181 


371 


676 


INC+SMM 








1.89 


184 


175 


282 


568 


with expansion 


8.1 


P 


2.1±0.2 


3.58 


183 


175 


488 


808 


INC+SMM 








2.0 


176 


158 


259 


529 


with expansion 


4.0 


4 He 


1.7±0.2 


3.89 


184 


177 


484 


836 


INC+SMM 








1.77 


177 


161 


238 


502 


with expansion 


14.6 


4 He 


2.2+0.2 


4.47 


173 


159 


723 


1132 


INC+SMM 








2.19 


154 


103 


183 


404 


with expansion 


22.4 


12 C 


2.2+0.3 


4.04 


163 


153 


924 


1216 


INC+SMM 








2.17 


139 


86 


207 


415 


with expansion 



Table 1 

The calculated properties of nuclear remnants from proj + Au collisions for 
INC+SMM and INC+Expansion+SMM. The M IM f is the mean number of IMF's 
for events with at least one IMF and Ar, Er are the mean mass number and exci- 
tation energy (in MeV), respectively, averaged over all inelastic collisions, while the 
quantities Amf, Emf are averaged only over residues decaying by IMF emission. 
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* p(2.1 GeV) 
o p(3.6GeV) 
▼ p(8.1 GeV) 
▲ a(4GeV) 

■ a(14.6GeV) 

• C(22.4GeV) 




100 120 
ENERGY(MeV) 



Fig. 1. Energy distribution of carbon fragments obtained for different collision sys- 
tems at = 89°. The lines are calculated in the INC+Expansion+SMM model 
assuming no flow. The spectra are arbitrarily scaled along the yield axis to avoid 
mixing of symbols. 
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Fig. 2. Upper part: The mean kinetic energies of outgoing fragments per nucleon 
measured at 6 = 89° for p(8.1 GeV), 4 He (14.6 GeV) and 12 C (22.4 GeV) collisions 
with Au. The shown error bars are due to statistics only, a systematic error of 5% 
has to be added. The lines are calculated using INC+Expansion+SMM approach 
assuming no flow. Lower part: Flow energy per nucleon (solid points) obtained as 
a difference of the measured fragment kinetic energies and the values calculated 
without flow. The dashed line shows the calculations (see text) assuming a linear 
radial profile of the expansion velocity with v® low = 0.1 c (at the surface). 
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<R Z >/R sys 



Fig. 3. Experimentally deduced mean flow velocities < Vfi ow > /c (left scale) for 
12 C+Au collisions (triangles) as a function of the fragment charge, and the mean 
relative radial coordinates of fragments < Rz > / R sys (right scale) , obtained under 
the assumption of a linear radial profile of the expansion velocity. 
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